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ABSTRACT 



Context. For the evolution of large-scale structures in the Universe, it is unclear whether active galaxies represent a phase that each 
galaxy undergoes, and whether and to which extent the evolution of black holes at their centers is important. Binary black hole (BBH) 
systems may play a key role in our understanding of the above questions. 

Aims. We investigate the Caltech-Jodrell Bank flat-spectrum (CJF) sample to identify evidence that supports the merger-driven evolu- 
tion scheme of active galaxies, and search for tracer-systems of active galactic nuclei (AGN) evolution and possible BBH candidates. 
We discuss the strength of and uncertainty in the evidence and formulate a set of selection criteria to detect such tracer-systems. 
Methods. We conduct an extensive literature search for all available multiwavelength data, particularly in the optical and infrared 
regime, and morphological information about the CJF sources. We perform a statistical analysis of the properties of this sample. 
Results. We find 1 ULIRG (Mrk 231) included in the CJF prototype of a transitory system. In total, 28.6% of the CJF sources with 
z < 0.4 are distorted or have a companion. Given the unbiased sample used here, this provides strong evidence of the ubiquity of 
the merger phenomenon in the context of active galaxies. We find a correlation between the radio and the near-infrared luminosity 
of the high-luminosity sources, interpreted in the context of the interplay between a star-formation and AGN component. We find a 
connection between variability and evolutionary transitory systems, as selected on the basis of their near-infrared colors. We select 28 
sources that trace the different evolution phases of an AGN, as well as a number of the most promising BBH candidates. We find 4 
sources with almost periodical variability in the optical and radio on similar timescales. 

Key words. Galaxies: statistics - Galaxies: active - Galaxies: nuclei - Galaxies: evolution - Galaxies: fundamental parameters 



1. Introduction 

As one of the most prominent predictions of the theory of 
General Relativity, black holes have been at the center of as- 
tronomical research for many years. There is a large amount 
of indirect evidence that most, if not all, galaxies host a su- 
perm assive black hole (SMBH) in their center (e.g.. lBlandfordl 
Il986l) . It is widely beheved that this SMBH is the main con- 
stituent of activity in active galactic nuclei (AGN) (for a review 
see e.g., Begelman et al. 1984). Although the mechanism that 
triggers this activit y remains unknown, it has been argued by 
many authors (e.g., Hutchings & Campbelliri983l:ISanders et al.l 



Il988t iHopkins et al.l 120061) that galactic mergers should play a 
predominant role in these processes. 

In the context of the hierarchical model (e.g., iFrenk et aP 

Il988l) . galaxies with a binary black ho le (BBH) at the i r cen- 
ter should be common syste ms (e.g., iMerritt & EkersI l2002t 
iHaehnelt & K aufFmannI |2002[). These systems, although pre- 
dicted in theory (Begelman et al. 1980; for a review on the evo- 
lution of BBH see Merritt & Milosavlievic 2005) have been dif- 
ficult to study. Only a few of them have been directly observed 
(e.g.. [Hudson et al.li2006t iKomossa et al.ll2003l; [Rodriguez et all 
l2006l) . while in a number of galaxies the pr esence of a BBH 
has been discussed (for an overview, see iKomossal l2006h . 
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Active galaxies hosting BBH will be primary candidates for 
the future gravitational wave space interferometer L.I.S.A. (e.g., 
Sesana et al. 2005). 

It is important to identify additional BBH candidates. In 
this paper, we use the Caltech Jodrell-Bank Flat-spectrum (CJF) 
flux-limited sample of 293 radio-loud active galaxies to search 
for these systems, while looking for evidence of the merger- 
driven evolution scheme for AGN and identifying AGN at dif- 
ferent evolutionary stages. This is achieved by means of a mul- 
tiwavelength analysis of the sample. 

This paper is organized as follows. In Sect.|2] we describe the 
CJF sample. In Sect. [3] we describe the different states of BBH 
systems, their properties, and the criteria to search for them in a 
large sample of sources. In Sect. 31 we present the criteria, trac- 
ers, with which we investigate the relevance of mergers in AGN. 
In Sect. |5] we present multiwavelength data for the CJF sample. 
In Sect. |6] we present the analysis of the data presented in Sect. 
|5] In Sect. 121 we present the candidates for different evolutionary 
stages and BBH systems. We conclude with a discussion in Sect. 
[HI a summary, and an outlook. Throughout the paper, we assume 
the cosmological parameters //q = 71 kms^^Mpc^^, Qm = 0.27, 
and Oa = 0.73. 

2. The CJF Sample 

The CJF sample (iTavlor et al.[[l996h consists of 293 radio-loud 



active galaxies select ed (see Table p from three different sam- 
ples (for details, see [Britzen et al.l[2007bl) . The sources span a 
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Table 1. CJF sample and its properties. 



Frequency(MHz) 


4850 


Spectral Index 


4850 > _o <; 
"1400 - "-J 


Declination 


5>35° 


Galactic latitude 


\b\ > 10° 


# Quasars 


198 


# BL Lac 


32 


# Radio Galaxies 


52 


# Unclassified 


11 


# Total 


293 



large redshift range (see Fig. 1 in iBritzen et al.ll2()08|). the far 



thest object being at a redshift z = 3.889 (1745+ 624:lHooket al 



19951) and the closest at z = 0.0108 (1 146+596: lde Vaucouleun 
199lh . The average redshift of the sample is Zavg = 1.254, 



ZBLLa, 



c,avg 



0.546, ZRCavg = 0.554, and zqso. 



a\'g 



1.489 for 

BL Lacs, radio galaxies, and QSOs, respectively. All the ob- 
jects have been observed with the VLBA and/or the global VLBI 
network. Each source has at least 3 epochs of observations and 
has been imaged and studied kinematica lly dBritzen et al.lll999t 
IBritzen et alir2007at lBritzenetalll2008h . The X-ray properties 
have been studied a nd correlated with their VLBI properties 
(IBritzen et al.ll2007bh . 



3. Binary black hole system properties and 
identification criteria 

Binary black hole systems are expected t o represent an inter - 
mediate evolutionary stage of, possibly all ( IBritzen et al.ll200ll) . 
active galaxies, that are the result of a prior merger (e.g., 
iMerritt & Milosavlievic 2005) . The evolutionary track of bound 
BBH systems has been studied extensively (for a review, see 
iMerritt & Milosavhevid l2005l) and is divided in three major 
phases: 

- as the two host galaxies merge, the two black holes sink to- 
wards the center of the merging system through dynamical 
friction; 

- the two black holes form a bound binary system and continue 
to lose energy by means of gravitational interaction with the 
central stellar population and its surrounding gas; 

- and finally, after reaching a critical separation, the systems 
rapidly lose energy by means of gravitational radiation, lead- 
ing to the coalescence of the two black holes. 

As bound (or "hard") BBH systems, we consider a system of two 
black holes for which its binding energy per unit mass exceeds 
the kinetic energy of the stellar population in the nucleus of the 
system. From the above, we expect the BBH properties to be di- 
rectly linked to the properties of the host galaxy and the activity 
of the resulting AGN. 

Given the distance of active galaxies and the timescales in- 
volved, one expects to observe BBH systems in the first two 
phases of their evolution, before or just after they have formed 
a bound pair. We are therefore interested in identifying and 
classifying systems that have recently or are curren tly undergo- 
ing a merger phase. According to simulations (e.g. jMaver et al.l 
12007 ). these systems exhibit complex morphologies because of 
the tidal forces exerted on the merging galaxies. 

Ultra luminous infrared galaxies (ULlRGs) probably host 
an early, unbound BBH system because they are most 
likely to be merging syst ems and the precursors of AGN 



3C 75. lKomossa et al.ll2003l:lHudson et al.l2006h . Almost half of 
the ULIRGs observed so far exhibit nuclear emission lines, typ- 
ical of LINER and Seyfert nuclei (e.g., Nagar et al. 2003]). All 
of th e above phenomena can be linked together (iHopkins et al.l 
l2006i) . illustrating the connection between galaxy, and AGN in 
particular, evolution and BBH research. 

Mergers are not the only systems in which searches for 
BBH systems have been performed. Binary black holes in 
the second stage of their evolution may reside in the nuclei 
of relaxed active galaxies, resembling to first order a single 
black hole AGN. Very Long Baseline Interferometry (VLBI) 
techniques have been used to map sources with complicated 
structures that have been explained i n terms of a BBH (e.g. , 
0420-014, Britzen et al. 2001; 3C345. Lobanov & Rolandll2005l: 
1803-H784. iRoland et al.ii2008i). A BBH system may pro duce a 
precessing jet base ("e.g.JYokosawa & In oue 1985; Romer oet al.l 
'2003'; 'RolandeLal.' 2008 ; Kaastr a &Rooslll992.) and thus a bent 
jet. The detected perio dicities in the hgh t curves may also be 
caused by a BBH (e.g.JLehto & Valtonen 1996; Sillanpa a et al.l 
Il988t lOian et al.1 120071: ICaproni & Abrahami i2004).' Finally, 
there is a number of more indirect ev idence possibly connected 
to these systems (e.g.. lKomossall2006h . 

4. Evolution criteria 

We search for evidence supporting the connection between ac- 
tivity in galaxies and merger events. We also identify and clas- 
sify candidate systems at different evolutionary stages using the 
following criteria: 

- merger signs (morphological distortions, tidal tails, compan- 
ions etc.); 

- excess infrared fluxes and starburst activity; 

- color information. 

In the context of the merger-driven evolution of AGN, we 
search in addition for BBH candidates . We assume that BBH 
systems are created in systems that have undergone a merger. 
It is expected that unrelaxed (morphologically disturbed) sys- 
tems host still unbound BBH systems, whereas bound pairs of 
black holes are expected to be foun d in more evolved sy stems 
(~ 10'' - 10^ years after the merger. iBegelman et allll980 ). We 
use variability as a tracer of these evolved systems. Although 
the above do not provide robust evidence of a BBH system, they 
may be used as strong filters for selecting possible candidates. 



4. 1 . Host galaxy signatures 

Merger signs - We search for morphological information that 
would characterize a system that is about to, is currently, or has 
recently undergone a merger event. We use this information to 
illustrate the connection between activity in galaxies and merger 
events. Systems with close companions are the precursors 
of merging systems and represent the earliest evolutionary 
stages of these systems. In this context, we collect in addition 
clustering information about our sources, as sources in clusters 
have a higher probability of taking part in galaxy interactions 
and mergers. Tidal tails and other morphological distortions are 
used as explicit signs of current or recent merger activity. 

Star-formation activity - A higher star-formation rate may 



^Ht 



dCanalizo & Stocktonll200Tl) . (as is the case of NGC 6240 and 



be caused by merger e vents (e.g..lBarne .s & Hernauist 1991) . It 
is known that ULIRGs (iSanders et aljl988|) are considered to be 
a direct result of merger events ("e.g. JSanders & Mirabe]|fl996l) 
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and are syste ms that are powere d by a strong starburst com- 
ponent (e.g., iGenzel et al.l 11998 1). although higher luminosity 
ULIRGs appear to be dominated by a nuclear component. We 
use starburst activity, as well as extreme infrared luminosities 
induced by intensive star formation, as additional evidence for 
systems possibly undergoing a merger. 

Color information - We use the available near-infrared (HJK 
photometry) information to identify systems that are possibly 
associated with merger events. Owing to their high concentra- 
tion of dust, merging systems appear redder in color than their 
non-merger counterparts or simple evolved stellar populations. 
We use this effect to identify evolutionary transitory systems. 

Infrared emission - Infrared emission from active galaxies 
originates partly in the host galaxy of the AGN, as a result of 
thermal emission processes (thermal radiation from heated dust 
grains and gas). Part of the infrared emission also stems from 
non-thermal processes, namely synchrotron emission, produced 
by the same seed population of electrons that produces the radio 
emission of the compact core and the jets of AGN. Infrared 
emission can be, partially, used as a tracer of the thermal 
emission of the AGN host galaxy. High resolution infrared 
observations can potentially resolve the AGN host and reveal 
any peculiar morphological features, signs of tidal interactions. 

Optical emission - The optical and UV continuum observed 
in AGN originate in the inner part of the galaxy, mostly from 
the accretion disk and its photosphere, and partly from the 
jet. The emission is predominantly thermal (as implied by the 
observation of the big blue bump. .Krolik..l999.) and, based on 
certain assumptions, the radiation resembles that of a black 
body, affected however by different absorption mechanisms. 
The non-thermal jet component may however dominate the total 
emission because of beaming effects. As in the infrared, any 
morphological disruptions may also be seen in the optical. 



4.2. AGN signatures 

Variability - Although variability can be induced by a number 
of different effects, almost periodical, multiwavelength vari- 
ations in the lightcurves of objects point towards a periodic 
modulation in the nucleus of the system. Assuming that a hard 
BBH system is possibly the origin of these variations, we use 
multiwavelength variability as a tracer of hard BBH systems 
and thus late-merger (or merger-remnant) systems. 

Radio Emission - Radio emission originates mostly in the 
compact core and jet of the active galaxy. It is produced by 
relativistic electrons, gyrating in a magnetic field, by means 
of synchrotron radiation. Radio properties of a source provide 
direct information about the properties of the AGN. The radio 
luminosity of an active galaxy can be used as a measure for its 
non-thermal, AGN component. 

A more detailed discussion of the above can be found 
later on in Sect. [8j 



5. Multiwavelength properties of the CJF sample 

We present multiwavelength information about the CJF sources 
gathered from our extensive search of the available literature. We 
focus on the radio, infrared, and optical parts of the spectrum. 



5. 1 . Radio emission 

The CJF is a flux-limited radio-selected sample of flat-spectrum 
radio-loud AGN. Although originally created to study the kine- 
matics of pc-scale jets and superluminal motion, the CJF repre- 
sents an ideal sample for an unbiased study of the evolution of 
active galaxies. Unlike previously used samples in the same con- 
text, the radio-selection of the CJF allows us to exclude biases 
towards transitory systems (as in the case of infrared-selected 
samples). 

The CJF sample (Ta ble[I]l has been most extensively studied 
in t he radio regime (e.g., Taylor et al. 1996; Pearson et al. 1998; 
Britzen et al.lll999l:j\%meulen et al..2003 : Poflacketal. 2003; 
Lowe et all 120071; iBritzen et all l2007bl: IBritzen et al. 2007a; 



Britzen et al.l l2008l) ~ lBritzen et alJ (l2008l) develop a localized 



method for calculating the bending of the jet associated with in- 
dividual components. The maximum of the distribution of local 
angles is at zero degrees, although a substantial fraction shows 
some bending (0 - 40 degrees). A few sources exhibit sharp 
bends of the order of > 50 degrees (see Fig. 13 in IBritzen et al.l 
2008). T hese morphologies on pc-scales (for a detailed discus- 
sion, see lBritzen et al.ll2008l) can be explained in the context of a 
BBH system. A future paper will study the morphologies of the 
jets of the CJF sources and their evolution. 

Although the CJF sample consists mostly of blazars, presum- 
ably highly beamed sources, the kinematical study of the sam- 
ple identifies a large number of sources with stationary, sublu- 
minal, or, at best, mildly superluminal outward velocities (e.g., 
see Fig. 15 in Britzen et al. 2008). Combined with a number 
of sources with inward ly moving com ponents (e.g., 0600-H422, 
1751-1-441, 1543+517. IBritzen et al]i2007 a). these sources do 
not fit into the regular paradigm of outward, superluminaly mov- 
ing components in blazar jets. One explanation of these pe- 
culiar ki nematic behaviors is that of a precessing, or helical, 
jet (e.g., iConwav & Murphvlll993h as a result of a BBH sys- 
tem. Other interpretations include traiUng shocks produced in 
the wake of a sing le perturbation propagating down the jet (e.g., 
lAgud o et al. 2001; but also see Mimica et al. 2009), or standing 
re-coUimation shocks (e.g.. Gomez et al. 1995) . 

Using the power spectral analysis method, iFan et al.] (l2007l) 
study the lightcurves of a large number of AGN at three dif- 
ferent frequencies, showing that there are different periodici- 
ties in them. There are 19 CJF sources in the Fan et al. sample 
(see Table 12]) that have almost periodical behavior in their radio 
lightcurves. Seven additional sources exhibit almost periodical 
variability in their radio lightcurves, as found in the literature (in 
total 26 sources, 8.9% of the sample; Tables |2] and [3] for refer- 
ences). Intra-day variability is not considered here (see Sect. [8] 
for details). 



5.2. Infrared emission 

There has been ongoing discussion about the evolutionary con 
nection between ULIRGs and AGN, with ULIRGs being as 
sumed to be the predecessors of active galaxies by means of a 
merging process (e.g., Canalizo & Stockton 2001; Bekki et aL 



2006; Ri echersetal.ll2008L iBennert et alJ 120081: iHopkins et al 



2006.) . We find that one so urce in ou r samp le can be classi - 
fied as an ULIRG (Mrk 231. ICao et a j. 2008; Kim et alj|2002h . 
There are 101 sources (34.5% of the sample; Table |4]l that have 
been observed and detected (we note that the presented data here 
are not statistically complete) in the ne ar to far infra red by the 
Infrared Astronomical SatelUte aRAS. ICleggil980l) and the 2 
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Table 2. CJF sources studied for variability across the spectrum. If available, we indicate the respective timescales. Columns (l)-(4) 
give the lAU source name, other name, the type of the source, and its redshift, respectively. Col. (5) gives the variability in the radio. 
Col. (6) in the infrared. Col. (7)in the optical. Col. (8) in the X-ray, and Col. (9) in the y-rays. For references, see Table[3] 



Source 



Other Names Type z 



Var^dio 

(yr) 



VariR 

(yr) 



Varopt 

(yr) 



Varxray 

(days) 



Var^ 

(days) 



0016+731 




Q 


1.781 


6.6 










0035+367 


4C 36.01 


Q 


0.366 






6.8'« 






0133+476 


OC457 


Q 


0.859 


5.8-10.9 










0212+735 




Q 


2.367 


15 










0219+428 


3C66A 


BL 


0.444 




yes 


4.25'°,2.5",6320'",65"-" 




39 
yes 


0316+413 


3C84 


G 


0.018 


9.7 




3.06, 65-68^2 


8003"* 




0402+379 


4C37.11 


G 


0.055 


yes^ 










0454+844 




BL 


0.112 






3 1232... 






0642+449 


OH+471 


Q 


3.396 






yes2* 


yes^^ 




0711+356 


OI318 


Q 


1.620 


0.7^ 




11.2^5 






0716+714 




BL 




5.4-15,5.5-6-* 


yes 10 


3.3'' 


147034,. 


4(1 

yes*" 


0804+499 


OJ508 


Q 


1.432 


1.1, 1.8,2.73 








yes3" 


0814+425 


OJ425 


BL 


0.245 


19.2-20.7 


yes" 


2,423,.. 






0836+710 


4C 71.07 


Q 


2.180 


5.7-9.0 








yes3'' 


0917+449 


VIPS 0214 


Q 


2.180 










yes3'' 


0923+392 


4C 39.25 


Q 


0.699 


22.5-22.4 


(?)'^ 


yes^^ 






0927+352 




BL 




6= 










0945+408 




Q 


1.252 


3.6-7.3" 










0954+556 


4C 55.17 


Q 


0.909 










yes3« 


0954+658 




BL 


0.368 






2623... 




yes3' 


1038+528 


OL564 


Q 


0.677 


9.3 










1101+384 


Mrk421 


BL 


0.031 


3.7-8.0 


yes».i4 


6.2'« 


6236 


39 
yes 


1254+571 


Mrk231 


G 


0.042 


yes' 




yes^" 


yes^' 




1418+546 


OQ530 


BL 


0.151 


2.5-11.3 


yes'3 


yes^' 






1504+377 




G 


0.671 




yes'5 








1624+416 


4C 41.32 


Q 


2.550 


16.7-37.1 










1633+382 


4C 38.41 


Q 


1.807 


14.3 




11 
yes 




yes3' 


1638+398 


NRA0 512 


Q 


1.660 






0.5-31^2 






1642+690 




Q 


0.751 


6.6-7.0 










1641+399 


3C345 


Q 


0.595 


5.^,3.5-4** 


yes'3 


5.6, 11.42*, iQi29^ 1930 


0.7234 




1652+398 


Mrk 501 


BL 


0.034 


6.8-12.3 


yes* 


yes27 


233* 


23.2"' 


1739+522 


4C 51.37 


Q 


1.381 


1.43 








yes3' 


1749+701 




BL 


0.769 


12.3 




yes^^ 






1803+784 




BL 


0.679 


7.0-11.3 




1123,.. 






1807+698 


3C 371.0 


BL 


0.051 




yes'" 


2.7'", 2.4'* 






1823+568 




BL 


0.664 


8.8-14.9 










1928+738 


4C 73.18 


Q 


0.302 


7.7-8.2 




yes22 






2007+777 




BL 


0.342 




yes" 


0.623", 10-403'" 






2200+420 


BLLac 


BL 


0.069 


3.9-7.8 


yes" 


7.832,0.6,0.88, 1433,7.7'* 




yes3'' 


2351+456 


4C 45.51 


Q 


1.986 


11.8 











Micron All Sky Survey (2MASS. lKleinmann|[T992l) . as well as 
other Earth-based telescopes. 

Only a few blazars have extensive light curves in the infrared. 
We identify 12 sources (4.1% of the sample) with infrared vari- 
abihty (Table [Hi. 

5.3. Optical emission 

Almost all CJF sources (285 sources, 97.3 % of the sample) have 
been detected in the optical (see Table 1 in iBritzen et al.ll2007bl 
for magnitudes and redshifts). 

We also investigate the CJF sample for the number of sources 
that belong to a cluster or a group of galaxies. Because of the 
large redshift span of the CJF, it is not possible to have cluster 
information f or sourc es at high z (z > 1.5, e.g., Brodwin et al. 
I2OO8 ; BlakeskeeLalj|2003). We find 54 sources (18.4% of the 
sample) belonging to a cluster or group of galaxies (see Table 



|6ll. For z < 1, we calculate that 25.0% of the QSOs (zavg = 
0.60), 33.3% of the radio galaxies {z„rg = 0.27), and 29.4% 
of the BL Lac objects {zavg - 0.25) are found in clusters. We 
note that radio galax i es app ear to be more often members of 
clusters. I Hill & Lillvl (Il99lh fi nd a similar trend in their red- 
shift range (e.g., see Fig. 1 1 of iHill & Lillvl[l99ll) . Twenty-two 
sources (7.5% of the sample) show variable fluxes in the optical 
(intra-day variability has not been taken into account). Some of 
these sources have extensive enough lightcurves to calculate the 
timescales of this variability (see Table|2]i. 



6. Analysis and results 

6.1. Mergers in fhe CJF 

Assuming that activity (whether nuclear or star-formation) in 
galaxies is indeed triggered by mergers, we expect to ob- 
serve morphological distortions in the optical. If these struc- 
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Table 3. References for Table |2l 



Symbol 


Reference 


Symbol 


Reference 


1 


Fan et al. (2007) 


22 


Bozvanetal. (1990) 


2 


Rodrisuez et al. f2006) 


23 


Heidt& Wagner (1996) 


3 


Kelly et al. (2003) 


24 


Pica & Smith (198J) 


4 


Raiteri et al. (2003) 


25 


Smith etal.(1993J 


5 


Machalski Sl Brandt (1996) 


26 


Kollatschnv et al. (1992) 


6 


Ciaramella et al. (2004) 


27 


Iwasawaetal. (1998) 


7 


Ulvestadetal. (1999a) 


28 


Webb etal. (1988) 


8 


Lobanov & Zensus (1999) 


29 


Zhang et al. (199^ 


9 


Abdo et al. (2009) 


30 


Babadzhanvants etal. (1995) 


10 


Haeen-Thorn et al. (2006) 


31 


Bvchkova et al. (2003) 


11 


Fan & Lin (1999) 


32 


Hagen-Thorn et al. (1997) 


12 


Neueebauer et al. (1986) 


33 


Fan etal. (1998) 


13 


Takaloetal. (1992) 


34 


Bian & Zhao (2003) 


14 


Gupta et al. (2004) 


35 


Zamorani et al. (1984) 


15 


Stickeletal. (1996) 


36 


Osone & Teshima (20(31) 


16 


Edelson & Malkan (1987) 


37 


Gallagher et al. (2002) 


17 


Peng et al. (2000) 


38 


Osone (2006) 


18 


Smith &Nair (1995) 


39 


Fan et al. (2002) 


19 


Belokon & Babadzhanvants (2003) 


40 


Chen et al. (2008) 


20 


Xie et al. (2002) 


41 


NeshDor (2000) 


21 


Lainela et al. (1999) 


** 


in days 



tures are observe d, the system must have recently (100-200 Myr, 
iLotz et al] l2008h undergone a merger. A typical merger life- 
time is of the or der of a few billion years. Using simulations, 
iLotz et al] ( l2008h . find that the peak in the morphological dis- 
turbances of the merging system happens just after the first pass 
of the two galaxies and after the actual merging of the nuclei of 
the merging galaxies. In-between, the system does not exhibit 
any quantitatively appreciable asymmetry. Merger-induced star 
formation can be used as evidence of merging systems in this 
intermediate phase. 

We find 16 sources (5.5% of the sample, 11 radio galax- 
ies, 4 BL Lac objects, and 1 quasar) exhibiting distortions 
in the optical and infrared (see Table [HJ. We find 14 sources 
(4.8% of the sample, 9 RG, 3 BL, and 2 Q) with one or 
more close companions (Table |8]l. The percentage of distorted 
sources or sources with comp anions is in good a greement with 
iHutchings & Campbelll ( Il983h and lSuracd d 19981) . who find that, 
for z < 0.4 a 30% of quasars show such properties (28.6% for the 
CJF). This underlines the ubiquity of these morphologically dis- 
torted systems and the importance of mergers in the AGN phe- 
no menon. The angle to the fi ne of sight assumed for radio galax- 
ies (lUrrv & Padovanilll995h allows the investigation of the host 
galaxy, unlike in the case of BL Lacs and quasars, where the jet 
dominates the emission. This may explain the higher percentage 
of radio galaxies in the distorted sources. Moreover, their lower 
average redshift compared to the quasar sample, contributes to 
the above effect. Finally, with lower nuclear luminosity, radio 
galaxy cores do not outshine their host galaxies, thus facilitating 
these morphological studies. 

Given the small number of CJF sources observed in the far- 
IR, information about the star-formation rates of these objects 
are scarce. Nevertheless, we fin d four sources tha t have signs 
of star-burst activity (0316+413, Richer et al. 1993; 1254+571, 
[Richards et al._2005.: 1652+398, 1807+6 98.Condon et al._200Z) 
and one additional candidate (0248+430 iDudik et alJl2005h . 

Finally, since we find one ULIRG (Mrk 231) in our sample, 
we can infer a percentage for the content of ULIRGs in the CJF 
sample. Counting only the IR detected sources, we have a 1.0% 
of the sample exhibiting very high luminosities in the far-lR (> 



10 Ljo/. This number is a lower limit, given that not all CJF 
sources have been observed in the far-IR or are observable at 
their respective redshifts. 

6.2. Luminosity correiations 

In Fig.[ri we plot the near-infrared (K band) luminosity against 
the K-corrected radio luminosity obtained by single dish obser- 
vations at 5 GHz. 

By measuring the infrared luminosity of our sources, we 
probe the stellar light of the active galaxies, as well as radia- 
tion from the dust and gas component of the sources. We are 
interested in the relation of this thermal component to the non- 
thermal one, given by the radio luminosity of the AGN. We find 
a broken correlation between these two quantities. The lower 
luminosity sources appear to have a near constant NIR lumi- 
nosity independent of the radio core luminosity. A standard 
linear regression infers a slope of 0.18 ± 0.12 with a correla- 
tion coefficient of 0.35. For higher luminosities, a linear depen- 
dence is found. A slope of 0.80 + 0.11 with a correlation coeffi- 
cien t of 0.79 is given by standard linear regression. Following 
Bro wne & Murphvl ( Il987h . we examine in addition the corre- 
lation between LradiolLNm and both Lradio and z. We find a 
stronger correlation with Lradio (correlation coefficient of 0.788) 
than to z (correlation coefficient of 0.498). This leads us to be- 
lieve that the correlation is not an artifact of luminosity evolu- 
tion. The large scatter of the data is expected, because of sig- 
nificant variability of active galaxies and the radio and infrared 
obs ervations not being contemporaneous. 

iRiekd (Il978h was first to confirm a correlation between NIR 
and radio continuum for a sample of Seyfert galaxies. This indi- 
cates that the infrared output of Seyfert galaxies is dominated by 
reradiation from dust in the host galaxy. In this context, we find 
that higher luminosity sources are more dominated by their AGN 
component (as traced by the radio luminosity), than lower lumi- 
nosity sources, for which an excess NIR luminosity is found. 
This excess is measured relative to the luminosity predicted by 
the linear dependency seen for higher luminosity sources and 
indicates that there is significant infrared output from star for- 
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Table 4. CJF sources observed and detected in the near infrared. Column (1) gives the lAU name of the source, Col. (2) the type of 
source (Q=Quasar, BL=BL Lac object, G=Radio Galaxy), Col. (3) the redshift z. Cols. (4)-(6) fluxes in three different NIR bands. 
Col. (7) the relevant references (Table |5]l, and Col. (8) gives the calculated spectral index in the NIR. 



Source 


Type 


z 


1.25yum 
(mJy) 


Near-IR 

1.65/jm 
(mJy) 


2.20yum 
(mJy) 


Reference 


onir 


0022+390 


Q 


1.946 






0.06 (0.03) 


42 




0035+413 


Q 


1.353 


0.64 (0.06) 


0.75 (0.13) 


1.37(0.24) 


43 


-1.29 


0108+388 


G 


0.669 






0.14(0.06) 


44 




0133+476 


Q 


0.859 


0.479 (0.02) 


0.813 (0.04) 


1.32(0.04) 


45 


-L74 


0145+386 


Q 


1.442 


0.40 (0.04) 


0.68(0.11) 


0.30(0.17) 


43 


0.55 


0153+744 


Q 


2.388 


0.93 (0.07) 


1.29(0.15) 


1.96(0.23) 


43 


-0.89 


0219+428 


BL 


0.444 


9.22 (0.25) 


12.5 (0.4) 


15.4(1.1) 


46, 11 


-0.83 


0249+383 


Q 


1.122 






0.12(0.03) 


42 




0251+393 


Q 


0.289 


0.70 (0.05) 


0.79(0.11) 


1.21 (0.19) 


43 


-0.89 


0316+413 


G 


0.018 


46.5 (0.7) 


57.8 (0.9) 


55.70 (0.8) 


47 


-0.10 


0344+405 


G 


0.039 


17.8 (0.6) 


26.0(1.0) 


20.0(1.0) 


48 


0.26 


0402+379 


G 


0.055 


19.1 (1.8) 


39.8 (2.6) 


40(3) 


43 


-0.38 


0546+726 


Q 


1.555 


0.84 (0.02) 


1.24(0.03) 


0.91 (0.02) 


43 


0.03 


0650+453 


Q 


0.933 


0.33 (0.02) 


0.49 (0.02) 


0.84 (0.02) 


43 


-1.58 


0651+410 


G 


0.022 


46.6(1.1) 


57.0(1.9) 


48.0(1.5) 


48 


0.03 


0707+476 


Q 


1.292 


13.2 (0.3) 


10.9 (0.5) 


74 (0.6) 


43 


LOS 


0710+439 


G 


0.518 






0.84(0.19) 


49 




0731+479 


Q 


0.782 


0.58 (0.02) 


0.52 (0.03) 


0.67 (0.02) 


43 


-0.17 


0733+597 


G 


0.041 


43.2(1.3) 


52.3 (1.9) 


45.1 (1.8) 


48 


-0.01 


0814+425 


BL 


0.245 


1.42(0.04) 


2.21 (0.06) 


3.14(0.07) 


11 


-L34 


0821+394 


Q 


1.216 


0.77 (0.02) 


1.15 (0.02) 


1.45 (0.02) 


43 


-1.07 


0824+355 


Q 


2.249 






0.10(0.06) 


50 




0831+557 


G 


0.240 


0.80 


1.03 


1.23 


51 


-0.72 


0833+416 


Q 


1.298 


0.36 (0.02) 


0.43 


0.35 (0.02) 


43 


0.09 


0836+710 


Q 


2.180 


1.00(0.02) 


1.30(0.03) 


1.50(0.02) 


43 


-0.68 


0847+379 


G 


0.407 






0.63 (0.04) 


50 




0923+392 


Q 


0.699 


1.22(0.06) 


1.20(0.12) 


1.66(0.22) 


43 


-0.53 


0927+352 


BL 








0.81 (0.09) 


50 




0945+664 


G 


0.850 


0.030 (0.001) 


0.005 (0.001) 


0.08 (0.01) 


15 


-1.89 


0954+658 


BL 


0.368 




3.17 




52 




1003+830 


G 


0.322 


0.86(0.13) 


2.0 (0.4) 


2.3 (0.4) 


48 


-L67 


1020+400 


Q 


1.254 


0.74 (0.04) 


1.00(0.07) 


1.53 (0.07) 


43 


-L27 


1030+415 


Q 


1.120 


0.07 


0.12 


0.30 


53 


-2.56 


1031+567 


G 


0.460 






0.19(0.13) 


49 




1038+528 


Q 


0.677 


0.65 (0.05) 


0.42(0.10) 


0.84(0.16) 


43 


-0.45 


1058+726 


Q 


1.460 


0.42 (0.02) 


0.44 (0.02) 


0.48 (0.03) 


43 


-0.20 


1101+384 


BL 


0.031 


71(5) 


81(8) 


83(15) 


11 


-0.25 


1124+455 


Q 


1.811 


0.57 (0.01) 


0.54 (0.02) 


0.66 (0.02) 


43 


-0.25 


1128+385 


Q 


1.733 






0.32 


53 




1144+402 


Q 


1.088 


1.17(0.06) 


1.69(0.15) 


2.53 (0.25) 


43 


-L35 


1144+352 


G 


0.063 


20.0 (0.8) 


34.0(1.3) 


23.3(1.5) 


44 


-0.24 


1146+596 


G 


0.011 


247 (2) 


311(5) 


248 (8) 


43 


0.03 


1226+373 


Q 


1.515 


0.33 (0.02) 


0.48 (0.02) 


0.41 (0.03) 


43 


-0.36 


1254+571 


Q 


0.042 


62.0(1.2) 


113(2) 


199 (3) 


54 


-2.04 


1306+360 


Q 


1.055 


0.23 (0.02) 


0.41 (0.02) 


0.54 (0.02) 


43 


-L49 


1309+555 


Q 


0.926 


0.65 (0.02) 


0.93 (0.03) 


1.36(0.03) 


43 


-1.28 


1333+459 


Q 


2.449 


0.34 (0.05) 


0.40(0.11) 


0.66 (0.20) 


43 


-1.15 


1347+539 


Q 


0.980 


0.47 (0.02) 


0.53 (0.03) 


0.65 (0.04) 


43 


-0.56 


1413+373 


Q 


2.360 


0.36 (0.02) 


0.36 (0.02) 


0.84 (0.03) 


43 


-1.50 


1418+546 


BL 


0.151 


6.68(0.31) 


9.2 (0.4) 


11.8(0.5) 


11,55 


-0.99 



mation regions that dominates over nuclear infrared emission. 
Therefore, higher redshift detected (eff'ectively, higher luminos- 
ity) sources are less likely to be found in a transitionary evolu- 
tionary stage, since their spectrum is exclusively dominated by 
their nucleus, effectively outshining and hiding any effects of an 
ongoing merger. 



6.3. AGN near-IR colors 

In the same context as above, we are interested in investigat- 
ing the interplay between the star-formation and nuclear compo- 
nent in the sources of our sample. For this reason, we calculate 
the near-infrared colors of the sources with available data, after 
corre cting all magnitudes (HJK bands) for Galactic extinction 
(from lSchlegel et al]|1998h . From these, we construct the near- 
IR color-color diagram seen in Fig.|2] The lines for a black body 



Table 4. Continued. 
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■nerger-driven 


evolution of active galaxies using 1 


the CJF sample 


Source 


Type 


z 




Near-IR 




Reference 


0!n1R 








l.25fim 


l.65fim 


2.20fim 














(mJy) 


(mJy) 


(mJy) 








1417+385 


Q 


1.832 


0.41 (0.02) 


0.52 (0.02) 


0.66 (0.02) 


43 




-0.83 


1432+422 





1.240 


0.02 




0.05 


56 






1435+638 


Q 


2.068 


1.10(0.02) 


1.00(0.03) 


1.40(0.03) 


43 




-0.42 


1442+637 


Q 


1.38 


0.51 (0.02) 


0.88 (0.02) 


0.71 (0.02) 


43 




-0.56 


1504+377 


G 


0.671 


0.06 (0.01) 


0.12(0.01) 


0.25 (0.02) 


15 




-2.62 


1531+722 





0.899 


1.25(0.12) 


2.0 (0.4) 


1.7(0.4) 


48 




-0.53 


1622+665 


G 


0.201 






3.3 (0.3) 


57 






1633+382 


Q 


1.807 






1.95 (0.07) 


49, 


58 




1637+574 


Q 


0.749 


1.70(0.02) 


1.90(0.02) 


2.50 (0.03) 


12 




-0.67 


1641+399 





0.595 


3.23(0.10) 


5.08 (0.25) 


7.5 (0.4) 


43 




-1.46 


1652+398 


BL 


0.034 


86.0 (2.0) 


115.0(2.0) 


99.0 (2.0) 


59, 


11 


-0.22 


1656+477 


Q 


1.622 


0.53 (0.02) 


0.60 (0.03) 


0.53 


43 




0.02 


1700+685 


G 


0.301 


0.64 (0.02) 


1.00(0.03) 


1.60(0.04) 


43 




-1.57 


1716+686 





0.339 


1.30 


1.70 


2.00 


43 




-0.73 


1719+357 


Q 


0.263 


0.83 (0.21) 


1.6(0.3) 


2.8 (0.4) 


48 




-2.10 


1726+455 


Q 


0.717 


0.61 (0.02) 


0.84 (0.03) 


1.01 (0.03) 


43 




-0.86 


1732+389 


Q 


0.970 


0.26 (0.02) 


0.49 (0.02) 


0.77 (0.02) 


43 




-1.87 


1744+557 


G 


0.031 


66(1) 


80(2) 


66(2) 


60 




0.05 


1749+701 


BL 


0.770 






2.2(1.7) 


55 






1755+578 


Q 


2.110 


0.33 (0.02) 


0.50 (0.02) 


0.67 (0.03) 


43 




-1.20 


1758+388 


Q 


2.092 


0.43 (0.02) 


0.38 (0.02) 


0.49 (0.02) 


43 




-0.21 


1803+784 


BL 


0.680 


5.50 (0.25) 


5.9 (0.4) 


7.1 (0.8) 


55,11 


-0.39 


1800+440 





0.663 


1.00 


1.20 


0.84 (0.02) 


43 




0.37 


1807+698 


BL 


0.051 


20.2 (0.9) 


25.4(1.1) 


29.7 (2.7) 


59, 


11 


-0.65 


1826+796 


G 


0.224 


1.10(0.20) 


1.3(0.3) 


1.9(0.4) 


48 




-0.89 


1823+568 


BL 


0.664 






4.98 


11, 


61 




1843+356 


G 


0.764 






0.13(0.09) 


44 






1946+708 


G 


0.101 


2.8 (0.3) 


2.4 (0.5) 


3.0 (0.5) 


57 




0.07 


2116+818 


Q 


0.084 


8.5 (0.4) 


12.5 (0.6) 


15.7 (0.7) 


48 




-0.95 


2200+420 


BL 


0.069 


38.4 (2.8) 


53(4) 


71 (14) 


43, 


11 


-0.80 


2351+456 





1.986 






0.63 


49 






2352+495 


G 


0.237 






0.9 (0.4) 


49 







Table 5. References for Table H 



Number 


Reference 


Number 


Reference 




42 


Sanchez & Gonzalez-Serrano (2003) 


52 


Wright etal. (1998) 




43 


Barkhouse & Hall (2001) 


53 


Carballo etal. (1998) 




44 


Snellen etal. (1996) 


54 


Sanders etal. (1988) 




45 


mooni et al. (1994) 


55 


Padovani et al. (2006) 




46 


Gonzalez-Perez et al. (2001) 


56 


Puschell et al. (1982) 




47 


Jarrett et al. (2003) 


57 


Villani & di Sereeo Alighieri 


,0999) 


48 


Skrutskie et al. (2003) 
de Vries et al. (1998) 


58 
59 


Xie etal. (1998) 




49 


Golombek etal. (1988) 




50 


Lillv etal. (1985) 


60 


KnaDD etal. (1990) 




51 


Heckman etal. (1983) 


61 


Kotilainen et al. (2005) 





and a power-law emitting body divide the color-color parameter 
space into three regions. Sources situated above the black-body 
line are predominantly thermal emitters. Sources beneath the 
power-law line are dominated by non-thermal emitting mech- 
anisms. For sources between the two lines, both components are 
important. 

In Fig. |2] we separated blazars from radio galaxies accord- 
ing to the different nature of these objects. We anticipate that 
radio galaxies are dominated by emission originating in their 
host galaxies. Indeed 8 out of the 1 1 radio galaxies depicted 
in Fig. |2] are situated above the black-body line. Radio galax- 
ies appear significantly redder in both J-H and H-K colors with 
(J - H)a,gjiG = 2.2+1.1 and {H - K)arg,Rc = 1.8 + 1.1. With 



the exception of nine sources, the remaining blazars are situated 
beneath that same line, most being situated around and beneath 
the power-law line. 

We are particularly interested in the sources found in the 
space between the two lines, hence showing transitory proper- 
ties from a thermal emission dominated state to a non-thermal 
dominated one. Indeed blazars in this transition area, with an av- 
erage H-K color of {H - K)avg - 0.8 14 + 0.027, are considerably 
redder than their equivalents found above the black-body line 
((// - K)avg - 0.19 + 0.05), sign of an excess of dust leading to 
the observed reddening. 

In the previous section, we inferred a possible connection be- 
tween multiwavelength variability, binary black holes, and hence 
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Table 6. CJF sources observed and detected in the optical regime as group or cluster members. Cols. (l)-(3) as in Table|4l Col. (4) 
gives the reference for clustering (Table |7]i. 



Source 


Type 


z 


Reference 


Source 


Type 


z 


Reference 


0010+405 


G 


0.255 


62 


1038+528 


Q 


0.677 


81 


0014+813 





3.366 


63 


1101+384 


BL 


0.031 


82 


0153+744 





2.388 


64 


1144+352 


G 


0.063 


83 


0212+735 





2.367 


64 


1146+596 


G 


0.011 


84 


0219+428 


BL 


0.444 


65 


1150+812 


Q 


1.250 


64 


0309+411 





0.134 


66 


1213+350 


Q 


0.857 


62 


0316+413 


G 


0.018 


67 


1216+487 





1.076 


64 


0636+680 





3.174 


68 


1254+571 


G 


0.422 


85 


0651+410 


G 


0.022 


69 


1415+463 


Q 


1.552 


76 


0700+470 


G 




70 


1435+638 


Q 


2.068 


65 


0710+439 


G 


0.518 


62 


1624+416 





2.550 


77 


0727+409 





2.500 


71 


1633+382 





1.807 


76 


0733+597 


G 


0.041 


72 


1636+473 





0.740 


76 


0821+394 





1.216 


73 


1637+574 


Q 


0.749 


86 


0831+557 


G 


0.240 


74 


1638+540 


Q 


1.977 


87 


0836+710 





2.180 


62 


1638+398 





1.660 


79 


0847+379 


G 


0.407 


75 


1642+690 





0.751 


64 


0850+581 





1.322 


73 


1641+399 


Q 


0.595 


88 


0859+470 





1.462 


76 


1652+398 


BL 


0.034 


89 


0923+392 





0.699 


62 


1716+686 


Q 


0.339 


64 


0945+408 





1.252 


65 


1744+557 


G 


0.031 


69 


0954+556 





0.909 


77 


1751+441 


Q 


0.871 


64 


0955+476 





1.873 


78 


1803+784 


BL 


0.680 


65 


1010+350 





1.414 


79 


1807+698 


BL 


0.051 


82 


1020+400 





1.254 


76 


2214+350 





0.510 


90 


1031+567 


G 


0.460 


80 











Table 7. References for Table |6l 



Symbol 


Reference 


Symbol 


Reference 


62 


Gelleretal. (1984) 


77 


Wu& Han (1995) 


63 


Chernomordik (1995) 


78 


Condon etal. (1990) 


64 


West (1991) 
_Dobrzvcki et al. (2002) 


79 
80 


Kim et al. (1991) 


65 


O'Dea etal. (1996) 


66 


Edge & Rottaerins (1995) 


81 


Hmtzen(1984) 


67 


Rylec&Windram(1968) 


82 


Nair(1997) 


68 


Cristiani et al. (1997) 


83 


White et al. (1999) 


69 


Zwickv et al. (1966) 


84 


Golev & Prugniel (1998) 


70 


LiDovka & Lioovka (2002) 


85 


Small et al. (2003) 


71 


Ledlow & Owen (1995) 


86 


Hintzen etal. (1983) 


72 


Zwickv etal. (1968) 


87 


Andernach et al. (1988) 


73 


T^son (1986) 


88 


Hutchings et al. (1998) 


74 


Dunn et al. (2005) 


89 


Miller et al. (2002) 


75 


Hill&Lillv(1991) 


90 


Gehren et al. (1984) 


76 


Maoz (J 995) 







post-merger or merger-remnant systems. We test this connection 
using the near-IR color properties of our sources. We divide our 
sources into three groups, as defined by the two lines seen in 
Fig. |2] We find that 42.9% of the sources situated between the 
black-body and power-law lines, effectively sources assumed in 
a transition from a merger (a system that is dominated by its ther- 
mal emission) to an active galaxy (a system that is dominated by 
its non-thermal emission), are variable at different wavelengths. 
Only 11.1% and 10.5% of the sources above the black-body line 
and below the power-law line, respectively, show variability. 



7. Tracing the evolution in the CJF 

One aim of this paper is to study the possible evolutionary stages 
of the AGN in the CJF sample. Different subsamples of the CJF, 
and the individual sources within them, can probe these different 
evolutionary stages. Mergers, starburst activity, and BBH sys- 
tems are integral parts of this study. 

One of the most prominent examples supporting the link 
between starburst and BBH systems is NGC 62 40, one of few 
sources directly observed to have a binary core dKomossa et alj 
l2003h . This system is an ULIRG, hosts two AGN, and 
also clearly exhibits ongoing starburst activity, making it an 
archetype of the evolutionary scenario discussed here. 

By selecting sources with a companion, one can investigate 
the earliest stages of the merging process. Accordingly, by se- 
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Table 8. CJF sources observed to have a companion and sources showing evidence of an ongoing interaction with another system. 
Columns (l)-(4) as previously, Col. (5) gives the optical magnitude. Col. (6) gives cluster information for comparison (Table |6]l. 
Cols. (7) and (8) give references in regard to companion galaxies and their separation respectively, and Col. (9) gives the references 
in regard to interacting systems. By a star, we indicate ambiguous evidence. 



Source 


Other Name 


Type 


z 


m 


Cluster 


Companion 


Separation 


Interaction 










(mag) 








(kpc) 




0018+729 




G 


0.821 












91 


0108+388 


OC314 


G 


0.669 


22 










91 


0316+413 


3C84 


G 


0.018 


11.9 


67 


92 




-600 


92 


0402+379 


4C 37.11 


G 


0.055 


18.5 




2 




29 


2 


0710+439 




G 


0.518 


20.4 


62 


93* 




53 


93 


0804+499 


OJ508 





1.432 


17.5 










97 


0831+557 


4C 55.16 


G 


0.240 


17.5 


74 


93* 




81 


93* 


0954+658 




BL 


0.368 


16.7 




94 




474 




1031+567 


COINS J1035+5628 


G 


0.460 


20.2 


91 


80 




65 


49 


1146+596 


NGC3894 


G 


0.011 


13 


84 








98 


1254+571 


Mrk231 


G 


0.042 


13.6 


85 


37 




361 


37 


1418+546 


OQ530 


BL 


0.151 


15.5 


82 








99 


1641+399 


3C345 





0.595 


16 


88 


22,95 




47 




1652+398 


Mkn 501 


BL 


0.034 


13.8 


89 








99 


1807+698 


3C 371.0 


BL 


0.051 


14.8 


82 


94 




81 


99 


1823+568 


4C 56.27 


BL 


0.664 


18.4 




61 




46 


61 


1843+356 


COINS J1845+3541 


G 


0.764 


21.9 




80* 




129 


91 


1946+708 




G 


0.101 


16.1 




96 




67 


96 


2021+614 


OW637 


Q 


0.227 


19 


80* 


80* 




36 




2352+495 


OZ488 


G 


0.237 


19 


80 


80* 




47 




91. Stanghellini et al. (1993) 


92. Colina & Perez-Foumon (1990) 


93. Hutchings et al. (198^) 






94.Stickeletal. (1993) 


95.Kirhakosetal. (1999) 


96. Perlman et al. 


(200 IJ 




97. Hutchinas et al. (1999) 


98. Pustilniketal. (2001) 


99. Pursimo et al. 


(2002) 
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Fig. 1. Luminosities in the near- infrared as a function of th e radio luminosity at 5 GHz, from single dish observations with the Green 
Bank radio telescope (data from'Gregorv & Co ndonI 19911) . Linear regression fits for the low radio luminosity (dashed line) and high 
radio luminosity (continuous line) CJF sources are shown. 



lecting the sources with disturbed morphologies and large in- 
frared fluxes, we probe an intermediate phase (e.g., Mrk 231). 
Finally, sources with relaxed morphologies but almost periodic 
variability at multiple wavelengths are likely to represent sys- 
tems in which an assumed BBH has sunk to the center of the 
system. 

iLotz et aP (l2008h conducted simulations of equal-mass gas- 
rich mergers, classifying the sources into six distinct merging 
phases: pre-merger, first pass, maximal separation, merger, post- 



merger, and remnant. Using this classification scheme for our 
sources, we classify objects into pre and post-merger stages. 

The incompleteness of the information available about our 
sample does not allow us to make unambiguous claims for the 
classification of some sources (especially when differentiating 
between for example either pre-merger and maximal separa- 
tion, or post-merger and remnant classes). Among the distorted 
sources, those with companions can be classified as either a pre- 
merger or maximal separation phase, while those with no com- 
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Fig. 2. Near-infrared color-color diagram for all CJF sources with available information. We plot the lines for a power-law emitting 
body (continuous line) and for a black body (dashed line). We note blazars with open symbols (circles for BL Lac objects and 
triangles for the quasars) and radio galaxies with filled ones (see text for details). 



panions can be categorized as being in either the first pass, the 
post-merger, or remnant phase. For the last group, this trans- 
lates, roughly, to an age of 2-4 Gyr after the initial approach of 
the progenitor systems. Exceptions to the above constitute four 
sources (3C 84, Mrk 23 1 , 3C 37 1 , and Mrk 501) that exhibit both 
distorted morpholog i es and considerable starburst activity (e.g., 
[Richards et"al]|2005l; ICondon et al.ll2002l) . These characteristics 
put them in the merger phase and would constrain their age to be 
< 2 Gyr after the initial approach. 

In Table |9] we indicate the 27 most suita ble representative s 
of diff'erent evolutionary phases (following Lotz et al. 2008J). 
Sources with companions but no sign of disturbed morpholo- 
gies, are selected as pre-merger systems. Sources with disturbed 
morphologies but no companions and no apparent starburst ac- 
tivity, are assigned to the first pass phase when the two systems 
are at a minimum separation for the first time. The next phase 
corresponds to that of the maximal separation between the merg- 
ing systems. For this phase, sources with disturbed morpholo- 
gies and detected companions are selected. Sources identified as 
undergoing the first pass may also be merger remnants. More 
information is required to make this distinction. 

In the context of BBH systems, the last two categories shown 
in Table|9]are the most interesting, as we expect the SMBH at the 
centers of the merging systems to be close enough for the grav- 
itational pull to bind them together. For the merger phase, we 
select sources that have disturbed morphologies and also show 
signs of ongoing starburst activity (similar to Mrk 231). For 
the final category we select sources that appear extremely vari- 
able across their spectrum (variable in at least 3 diff'erent wave- 
length regimes) in an almost periodic manner. These sources 
do not show signs of ongoing interaction, merging (apart from 
1418-1-546, see Individual Sources), or starburst activity. These 
are probably BBH systems. In the context of the above catego- 



rization, we can identify these sources as possible post-merger 
systems. 

At least three of the sources in Table |9] (0954-1-658, 
0108+388, and 0248-H430) cannot be unambiguously catego- 
rized as one evolutionary phase, as they exhibit properties com- 
mon to more than one stages of evolution. Sources selected in 
the last two categories (merger and post-merger) are the most 
likely BBH candidates in the CJF sample. We plot their NIR col- 
ors in Fig. [3] As expected, most candidate sources are found to 
be above the power-law line, while most post-merger classified 
sources (selected solely by their variability properties) are found 
just above the power-law line, in the transitory area of the color- 
color diagram. This provides us with an additional selection tool 
for transitory systems (see Sect.lAt. 



8. Discussion 

In previous sections, we have presented the accumulated litera- 
ture information available for the CJF sources and a sample of 
CJF sources that trace AGN evolution, indicating the most likely 
BBH candidates among them. We now investigate the criteria 
used to identify candidates as BBH systems. It is interesting to 
discuss which subset of these criteria can provide us with a ro- 
bust enough argument for the existence of a BBH. 



8.1. Anticipated AGN properties 

Before we proceed with the discussion, we are interested in as- 
sessing the completeness of the information for the CJF sample 
presented here, as these were collected from the available litera- 
ture and are not the result of a uniform observational campaign. 
To achieve this goal, we create two subsamples of CJF sources. 
For the first one, we selected sources with available near-infrared 
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Table 9. CJF sources categorized in different evolutionary phases of AGN (foUowing lLotz et ani2008L for details see text). Columns 
(l)-(4) as in previous tables, Col. (5) provides information about companions. Col. (6) gives morphology information. Col. (7) 
indicates the number of wavelength regimes the source has been observed and detected as variable, and Col. (8) describes the 
starburst activity. For references, see previous tables. 



Source 


Other Name 


Type 


z 


Companion 


Morphology 


Variability 


Starburst 


Pre - merger 


0954+658 




BLLac 


0.368 


yes 




2 




2021+614 


OW637 


Q 


0.227 


yes 








2352+495 


OZ488 


G 


0.237 


yes 








First Pass 


0018+729 




G 


0.821 




disturbed 






0108+388 


OC314 


G 


0.669 




disturbed 






0804+499 


OJ508 


Q 


1.432 




tidal tail 


1 




1146+596 


NGC3894 


G/LINER 


0.011 




interacting 






Max Separation 


0402+379 


4C 37.11 


G/Sy 


0.055 


yes 


mteractmg 


1 




0710+439 


01+417 


G/Sy 


0.518 


yes 


disturbed 






0831+557 


4C +55.16 


G/LINER 


0.240 


yes 


disturbed 






1031+567 




G/Sy2 


0.460 


yes 


disturbed 






1823+568 


4C 56.27 


BLLac 


0.664 


yes 


non-relaxed isophotes 






1843+356 




G 


0.764 


yes 


disturbed 






1946+708 




G 


0.101 


yes 


non-relaxed isophotes 






Merger 


0248+430 




Q 


1.310 








possible 


0316+413 


3C84 


G/Sy2 


0.018 


yes 


twisted isophotes 


3 


yes 


1254+571 


Mrk231 


Q 


0.042 


yes 


disturbed 


3 


yes 


1652+398 


Mrk501 


BLLac 


0.034 




extra nuclear component 


5 


yes 


1807+698 


3C 371.0 


BLLac 


0.051 


yes 


interacting 


2 


yes 


BBH/Post 


•merger 














0219+428 


3C066A 


BLLac 


0.444 






3 




0716+714 




BLLac 








5 




0814+425 


OJ425 


BLLac 


0.245 






3 




0923+392 


4C 39.25 


Q 


0.699 






3 




1101+384 


Mrk421 


BLLac 


0.031 






5 




1418+546 


OQ530 


BLLac 


0.151 




extra nuclear component 


3 




1633+382 


4C 38.41 


Q 


1.807 






3 




1641+399 


3C345 


Q 


0.595 


yes 




4 




2200+420 


BLLac 


BLLac 


0.069 






4 





information, as this information was primarily presented and an- 
alyzed in this paper In addition, we excluded sources without 
redshift information. We identified 74 of these sources, the num- 
ber of sources being explicitly controlled by the availability of 
infrared information (see Table H)). For the second subsample, 
we selected those CJF sources that have been observed and have 
available information across the whole spectrum (radio to y- 
rays). We found only 10 of these sources, the number of sources 
mainly being regulated by availability of y-ray information (see 
Tablelinii. 



8.1.1. Near-IR subsample 

This subsample, having been selected on the basis of detec- 
tion in the near-IR, is expected to be more sensitive to in- 
frared bright sources and consequently have a higher probability 
of containing sources associated with merger events. This bias 
should be accounted for We begin by investigating whether this 
subsample is differs fundamentally from the CJF sample from 
which it was selected. The subsample has an average redshift of 
Zavg,sub - 0.864, effectively probing more nearby sources than 
the CJF (Zavg,cjF = 1.254). We also check the radio spectral in- 
dex af/AA, a measure of the compactness of the sources. The sub- 



Table 10. The CJF source with available information across the 
electromagnetic spectrum (radio, near-infrared, optical. X-ray, 
and y-ray). 



Source 


Other Name 


Type 


Redshift 


0133+476 


DC 457 


HPQ 


0.859 


0219+428 


3C066A 


BLLac 


0.444 


0316+413 


3C84 


G 


0.018 


0650+453 




Q 


0.933 


0814+425 


OJ425 


BLLac 


0.53 


0836+710 


4C 71.07 


LPQ 


2.18 


1101+384 


Mrk421 


BLLac 


0.0308 


1633+382 


4C 38.41 


LPQ 


1.807 


1652+398 


Mrk 501 


BLLac 


0.03366 


1803+784 




BLLac 


0.6797 


2200+420 


BLLac 


BLLac 


0.0686 



sample has lower average and median values of the index than 
the whole sample. This indicates that this subsample contains 
less compact objects. Finally we check the range of luminosi- 
ties for which almost complete information for the CJF is avail- 
able (in radio, optical, and X-ray) for both the subsample and the 
CJF as a whole. In both samples, we find that the ranges are the 
same for all three different wavelengths. We conclude that this 
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Fig. 3. Near-infrared color-color diagram for the sources in- 
cluded in Table |9] with available information. We plot the lines 
for a power-law emitting body (continuous line) and for a black 
body (dashed Hne). We note two objects (Mrk 231, NGC 3894) 
of interest due to their transitionary nature (see text for details). 



subsample probably contains sources that are closer and more 
extended but is otherwise quite similar to the CJF. 

In the context of merger-driven evolution, we find that 18.9% 
of the sources in this subsample, exhibit morphological distor- 
tions and starburst activity. Eleven sources (14.9%) are found to 
have companions. Compared to the CJF sample in total (5.5% 
of sources exhibit distortions; 4.8% have companions), merger 
effects appear to be more prominent in this subsample, as ex- 
pected. However, the more than three times higher percentage of 
distorted sources and sources with companions than the CJF as 
a whole indicates that the true number of these sources might be 
higher than currently believed. 

We are also interested in the variability properties of this sub- 
sample: 22.7% of the sources are variable in the radio, 13.3% in 
the infrared, 18.7% in the optical, 6.7% in the X-rays, and 12.0% 
in the y-rays. For most of these sources and wavelength regimes 
(with the exception of infrared and y-rays), periodicities have 
been found in their lightcurves (see Tables|2]and[3]for timescales 
and references). Compared to the total CJF sample, these sources 
appear to be significantly more variable. This is an unexpected 
result, as nothing in our selection criteria and the properties of 
this sub-sample implies a bias towards more variable sources. 

8.1.2. Multiwavelength subsample 

This second subsample was selected to assess how our results 
(concerning the importance and relevance of merging events and 
variability for active galaxies) is affected by the availability of 
multiwavelength information (or lack thereof). Given the very 
small number of the sources included in this subsample, a sta- 
tistical comparison to the whole sample is difficult. The average 
redshift of this subsample is Zm-g.sub - 0.689, and the subsample 
consists almost entirely of blazars (with the exception of 3C 84). 
In the context of merger-driven evolution, we find that for 
this subsample, 17.6% of the sources (all of which are at z<0.05) 
exhibit morphological distortions and starburst activity. One 
source (9.1%) is found to have companions. Compared to the 
CJF sample in total, merger effects appear to be more promi- 
nent in this sub-sample. One must again take into account that 



the average redshift of the CJF sample is Zavg.cJF - 1 -254, much 
higher than this subsample. However, the more than three times 
higher percentage of distorted sources compared to the CJF as a 
whole indicates that the true number of these sources is higher 
than currently believed. 

The variability properties of this subsample are: 81.8% of the 
sources are variable in the radio, 45.4% in the infrared, 72.7% 
in the optical, 27.3% in the X-rays, and 81.8% in the y-rays 
(see Tables |2] and [3] for timescales and references, respectively). 
Compared to the total CJF sample, these sources appear to be 
highly variable, in an almost periodical manner This can be ex- 
plained as a bias effect, as sources that are highly variable in 
one waveband are more likely to be targeted for observations in 
other wavelengths. Alternatively, one can argue that (periodical) 
variability is an ubiquitous phenomenon for a large portion of 
active galaxies. The latter is also supported by the high percent- 
ages of variable sources detected for the first sub-sample. Multi- 
wavelength variability, however, can only be detected by exten- 
sive (and therefore time-costly) monitoring of a source, thus ef- 
fectively reducing the number of known sources with such prop- 
erties. 

For both subsamples, we find rather high percentages of 
sources with signs of past, present, or future merger events, as 
well as sources with variable fluxes at different wavelengths. 
Although the criteria used to select these two subsamples were 
different, the resulting percentage of sources with merger signs 
are the same for both subsamples, significantly higher than the 
percentage for the whole CJF. This implies that because of the 
incompleteness of the information available to us (due to both a 
lack of observations and instrumental limitations), the percent- 
age of active galaxies found to be linked to a merger event in our 
sample is only a portion of the true number, thus is only a lower 
limit. If the percentage calculated for the two subsamples were 
closer to the true percentage, we would expect approximately 50 
CJF sources to have been found taking part in a merger event 
and 31 sources having close companions. From the available in- 
formation, we have found 15 and 1 1 sources, respectively. 

8.2. Periodicities in fhe iightcurves 

Having concluded that this statistical analysis of the CJF can 
provide us with only lower limits for, e.g., the percentage of 
mergers in our sample of active galaxies, we can now discuss 
the criteria used to identify BBH systems. 

Periodicities in the lightcurves can be linked to a pe- 
riodic process at the center of an AGN, presumably a bi- 
nary system with a quasi-stable orbit that produces periodi- 
cal p erturbations in the accretion disk and the jet of the AGN 
^(e.g., iLobanov & Roland 2005: Oian et al. 2007.; Riegei] 120041; 
Roland et al. 2008; Britzen et al. 2009, in press). However, this 
is not the only possibility. 

Precession of the jet has been employed to explain a variety 
of phenomena includi ng the variability that AGN sh ow at differ- 
ent w avelengths (e.g.. I Caproni & Abraham.2004: .Roland et al.l 
Il994l) . An intrinsically helical jet can also exhibit the same mor- 
phological and kinematical results (for a study of the different 
kind of helicities in AGN jets, see Steffen 1997). Helicity in 
AGN jets may be the result of internal Kelvin-Helmholtz insta- 
bilities and in particular the m - 2 helical mode of the K-H in- 



stabilities (for a 2-D tre atment of the problem see e.g..|Batv et al 
2003t alternatively, see lPerucho et alJ[2006HMeier & Nakamura 
200i. 

Doppler beaming may also cause the observed variability in 
the lightcurves of AGN. Camenzind & Krockenberger (1992) 
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first described the so-called lighthouse effect. This model has 
been successfully applied to different sources (e.g., 3C 273, 
Came nzind & Krockenberger 1992; 3C 345, ISchramm et alj 



T993h . 



Different physical phenomena are described by different 
characteristic times. The BBH and precessing jets can de- 
scribe periodicities of the order of 10 years (see references 
above), whereas lighthouse models tend to explain the intra- 
day variability, also observ ed in many blazars (for a review, 
see I Wagner & Witzell 119951) . Intra-day variability can also be 
explained by interste llar scintillatio n (for a review, see e.g., 
iRick ett 1990; also see ' Melros3l 19941) . in addition to BBH mod- 
els (iRoland et al. 2009). 

Correlated variability at different wavelengths implies that 
the same population of particles produces this variability, thus 
pointing to an intrinsic effect, rather than, for example, a tur- 
bulent screen between the observer and the observed source. 
One of the most prominent sources displaying correlated vari- 
abilit y at different wavelengths is 0716-H714 (rad io and opti- 
cal) to uirrenbach et al.lll99ll; IWagner et alj|1996h . The source 
OJ287 is another famous example, that has c orrelated v ariabil- 
ity in the radio, optical, and infrared (see e.g.. Fan & Lin 1999^ 
OJ28 7 is a prominent candidate BBH system ( Sillanp aa et al 
19881: iLehto&Valtonen 1996; Pietila 1998; Valton en et al 



19991; iLiu & Wuii2002tiValtonenii2007i) . 



We categorized possible post-merger candidates as sources 
that exhibit almost periodic variabilities in their lightcurves 
(Table |9l), assuming that a BBH causes the variability at differ- 
ent wavelengths. Indeed we find that variability is connected to 
transitory systems as they are defined on the near-lR color-color 
diagram. In Table|2] we indicate the variability timescales where 
available. Out of 9 post-merger classified sources, 4 sources 
(0716-1-714, llOl-t-384, 1641+399,220+420) show almost peri- 
odic variability with similar timescales (~ 10 years) in the radio 
and optical. These are candidate hard BBH systems. More de- 
tailed investigation of their radio jet morphology may provide us 
with additional information. 



8.3. Infrared emission, merging scenarios, and starbursts 

The infrared emission from an AGN consists of different compo- 
nents, coming from different parts of the nucleus. Distinguishing 
and separating these constituents (accretion disk, thermal, and 
jet, non thermal) is not straightforward. Both components would 
be affected by a BBH system, since both would precess under 
the gravitational influence of the binary, although possibly on 
different timescales. 

In light of the merger-driven evolutionary scenario, we in- 
vestigate whether there is a connection between an excess in 
the infrared regime and sources with companions or sources 
that exhibit disrupted morphologies. For 16 sources with pecu- 
liar morphologies, we calculate an average infrared (K band) 
luminosity log{vLy)NiR^avg - 44.17 + 0.14 ergjs for an aver- 
age redshift of Zavg - 0.41. We also calculate the average lu- 
minosity in the NIR for sources that are found to have com- 
panions log{vLy)NiR,avg - 44.61 + 0.21 erg Is, for 11 sources 
with a Zavg - 0.34. Defining a subsample of the CJF sources 
with similar average and median redshifts (zavg - 0.38), we find 
log{vLy)MiR,avg - 44.5 + 0.1 ergjs. Sources with disturbed mor- 
phologies are of below average near-infrared luminosity, in con- 
trast to sources with companions that are on average and within 
the eiTors, equally luminous to the rest of the sample. The avail- 
able data do not show a prominent near-lR excess for AGN with 
perturbed morphologies as would be expected. The lack of mid- 



and far-infrared data does not allow us to extend this over the 
whole infrared range. 

Mrk 231 is an ULIRG, classified as a broad line quasar 
(BLQ) and a type 1 Seyfert. The system shows signs of an 
ongoing interaction and merger (ICanahzo & Stockton! 120011: 
iPetrosian et al.ll2002h . The galaxy hosts an AGN, with a Seyfert 
1 spectrum, but is however also classified as a quasar. The 
source exhibits strong starburst activity (R ichards et al . 2 005i) , 
which indicates that it has undergone a merger. The source has 
been studied in X-rays, where it exhibits extended emission, 
in addition to a hard spectrum central source that remains un- 
resolved and is off-center to the extended X-ray halo. An off- 
nuclear source has be en detected at X-rays at a distance of 6" 
(Gallag her et al.lf2002h . Mrk 231 has been studied for variability 
at different wavelengths. It has been found to be variable in the 
X-rays ([G allagher et al. 2002), optic al (in the ab sorption lines; 
lKollatsclinvetal.1 11992). and radio (lUlvestad etal.l l 1999b). 

Mrk 231 was observed in three epochs (Sep 1992, Sep 
1994, and Aug 1996) with the VLBA. The so urce remained 
unresol ved during all three epochs (see Fig. |4]l. lUlvestad et al.l 
(Il999bh and|Ulvestad et al. ( 1999a) however detected a pc-scale 
jet consisting of on e component moving at sub-luminal speeds. 
ICariUi et aD(ll998l) present strong arguments for the existence of 
a sub-kpc molecular disk formed in Mrk 23 1 . The above proper- 
ties point to a young AGN, created by a recent merger that has 
formed a molecular disk fueling gas into the nucleus of the newly 
formed AGN and a young radio jet, whic h agree well with th e 
merger-driven evolution scenario (see also. lHickox et al.ll2009l) . 

We note that 1 146+596 is somewhat similar to Mrk 23 1 . This 
low redshift (z=0.01 1) radio galaxy appears to be extremely lu- 
minous in the near infrared ( Skruts kie et al. 2006) . It shows signs 
of ongoing interaction (Pustilni ket al.ll2001[) to which t he in- 
frared excess can be attributed. Moreover. lPeck & Tavlorl (Il998h 
detect redshifted H I gas in absorption, which appears to repre- 
sent infall of gas towards the central engine. 

In the above context, for the merger phase (in Table |9]l we 
select sources that appear disturbed and also show signs of ongo- 
ing starburst activity (Mrk 231, 3C 84, Mrk 501, 3C 371). These 
are primary candidates for BBH that have not yet sunk into the 
center of the merging system (similar to NGC 6240). Most of 
these sources exhibit variability at different wavelengths. Three 
sources (0316+413, 1418+546, 1652+398) show almost peri- 
odic variability in the radio. 0316+413 and 1652+398 are in ad- 
dition members of clusters and therefore good candidates to be 
multiple merger systems. 

We note here the ambiguity of the classification of sources 
into the different pre-merger stages (Table |9j. Starburst activity 
is the defining, although not unambiguous, property of identify- 
ing systems that are currently undergoing a merger, the alterna- 
tive evidence discussed above being a weaker selection criteria. 



8.4. Binary black holes in the CJF 

We find 9 sources in the CJF sample that have been previously 
argued to be good BBH candidates. We present these sources in 
Tables nn and [TJl along with the relevant references. 

We compare the luminosities of these systems individually 
with the average luminosities at different wavelengths of a sub- 
sample of the CJF with similar median redshifts. All 9 sources 
are found to be significantly more luminous in X-rays. Two of 
these 9 sources (4C 37. 1 1, BL Lac) are found to have above av- 
erage K-band luminosity, possibly indicating a younger age for 
these sources. These two sources are fainter in X-rays than the 
remaining 7 sources. According to lSanders et alJ (Il988h . after a 
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Fig. 4. Three epochs of VLBI ob servations at 5 GHz o f the object 1254+571, Mrk 231. Contour maps with maximum flux at 0.143 
Jy/beam. Maps reproduced from lBritzen et alJ (l2007ah . 

Table 11. Kinematic and morphological information for CJF sources already considered to be BBH candidates. Columns (l)-(4) 
as in previous tables. Col. (5) shows the sources with X-ray detected jet. Col. ( 6) gives the misalign ment angle between the pc and 
the kpc-scale jet of the source. Col. (7) indicates the bending as calculated in Britze n et al.l (l2008h for the inner-most component 
identified in the jet. Col. (8) gives references related to the environment of the source such as a companion or an interaction with 
other systems. Col. (9) comments on the general kinematics of thepc-scalejet. Col. (10) gives the reference for possible precession 
models for the source, and Col. (11) gives the references in regard to possible BBH. For references, see Table [12] 



Source 


Other Name 


Type 


z 


X-ray Jet 


APA 

(°) 


Bending 

(o) 


Environment 


Kinematics 


Precession 


BBH 


0316+413 


3C84 


G 


0.018 




20 


51.9 


92 


inward"" 




108 


0402+379 


4C 37.11 


G 


0.055 




27/12 




2 


stationary '"- 




2,102 


0716+714 




BL 


0.300 




75 


2.6 




stationary '"^^ 


106 


109 


1101+384 


Mrk 421 


BL 


0.031 




13 


1.9 


82 


slow'""* 




110 


1641+399 


3C345 





0.595 


knot 


98 


46.8 


88,95 




107 


111 


1652+398 


Mrk 501 


BL 


0.034 




83 


10.5 


89,99 


stationary'"'* 


106 


110 


1803+784 




BL 


0.680 




70 


18.0 


99 


stationary'"'' 




112 


1928+738 


4C 73.18 


Q 


0.302 


knot 


22 


23.4 




stationary"" 


106 


113 


2200+420 


BLLac 


BL 


0.069 




30 


24.8 


100 






114 



merger has taken place, the peak of the SED of the system is 
pushed towards higher frequencies, as the AGN blows away any 
obscuring material remaining from the merging phase. This is 
another indication that 4C 37.11 and BL Lac are indeed in an 
earlier evolutionary stage than the other BBH candidates. 

In addition, we note that all BBH candidates exhibit period- 
icities across the spectrum. We find that 7 out of the total of 9 
sources have bent VLBA jets. Out of the sources in Table (TT] 6 
sources belong to a cluster or have a companion, illustrating the 
higher probability of these sources to harbor a BBH system. 



9. Conclusions 

The focus of this work has been twofold: (i) to investigate the 
evolutionary connection between active galaxies and merger 
events, and (ii) identify AGN candidates at different evolutionary 
stages, including binary black hole system candidates. To this 
end, we collect and present multiwavelength information about 
the CJF sources. 

From the previous analysis, it becomes apparent that merg- 
ers or merger-induced effects (starburst activity) are ubiquitous 
in active galaxies. With almost 30% of our sample showing ev- 
idence of interactions at z< 0.4 and with indications that more 
such systems exist, although not yet observed, mergers appear 
to play a predominant role in the evolution not only of galaxies 
but AGN as well. We note here the importance of the absence 



of any biases towards transitory systems in the selection crite- 
ria of our sample. This supports a direct link between merger 
events and active galaxies in general. Furthermore, we show that 
multiwavelength variability may be linked to transitory systems. 
Assuming that multiwavelength variability originates in binary 
black hole systems, this provides an additional link between 
merger events, binary black holes, and active galaxies, as well 
as a possible identifier of these sources. 

Following the merger-driven evolution scheme, we have also 
presented a first subsample of 28 CJF sources that are good trac- 
ers of evolution in AGN. We find at least 14 promising BBH can- 
didates at different points in their evolution (5 in an unbound and 
9 in a hard state). Classification of these sources cannot be un- 
ambiguous, given the lack of complete information. Additional 
observations are needed to clarify their true nature. 

Following arguments from both analytical and numerical 
simulations, it is anticipated that binary black holes should form 
during merger processes. We show that at least 13.6% of the CJF 
sources exhibit variability at different wavelengths (not includ- 
ing IDV), most of them being almost periodic. We show that 
most of the previous BBH candidates are highly variable across 
the spectrum, showing the tight connection between the two phe- 
nomena. The morphology and evolution properties of the AGN 
jets can also be linked to the evolution scheme and BBH systems. 
The morphology of CJF jets and its evolution will be studied in 
detail in a future paper. 
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Symbol 


Reference 


Symbol 


Reference 


100 


Disnev (1974) 


108 


Pronik(1988) 


101 


Britzen et al. (2007a) 


109 
110 


Nesci et al. (2005) 


102 


Maness et al. (2004) 


De Paolis et al. (2003) 


103 


Britzen et al. (2006) 


111 


Lobanov & Roland (2005) 


104 


Kellermann et al. (2003) 


112 


Roland et al. (2008) 


105 


Britzen et al. (2008) 


113 


Roosetal.(1993) 


106 


Lu& Zhou (2005) 


114 


Zier (2007) 


107 


Caoroni & Abraham (2004) 







Table 13. Variability information and, if available, timescales concerning CJF sources already considered to be BBH candidates. 
Col. (4) in the radio. Col. (5) in the infrared. Col. (6) in the optical. Col. (7) in the X-rays, and Col. (8) in the y-rays. Columns 
(l)-(3) as before. For references, see previous tables. 



Source 


Type 


z 


VarRadio 

(yr) 


VariR 


Varopt 

(yr) 






Varx 

(days) 


Var^ 

(days) 


0316+413 


G 


0.018 


9.7 














0402+379 


G 


0.055 


yes 














0716+714 


BL 


0.3 


5.4-15 


yes 


3.3 








yes 


1101+384 


BL 


0.0308 


3.7-8.0 


yes 


6.2 






62 




1641+399 


Q 


0.595 


5.1,3.5-4 


yes 


5.6, 11.4, 


10.1 


, 1.9 






1652+398 


BL 


0.03366 


6.8-12.3 


yes 


yes 






23 


23.2 


1803+784 


BL 


0.6797 


7.0-11.3 




0.003 










1928+738 


Q 


0.302 


7.7-8.2 




yes 










2200+420 


BL 


0.0686 


3.9-7.8 


yes 


7.8, 0.6, 0.88, 


14, 7.7 







More specifically, we find that: 

- one CJF source (Mrk 231) has been classified as an ULIRG. 
Its properties fit well in the merging paradigm of AGN. We 
also note that the similar source 1 146+596, which has a high 
luminosity in the near-infrared, is probably a merging sys- 
tem. Both sources are primary BBH candidates. 

- a coiTelation between radio luminosity and NIR luminosity 
exists. This is qualitatively explained in terms of reradiation 
of the non-thermal AGN emission. Low luminosity AGN ap- 
pear to have a near-constant NIR luminosity. 

- 4.4% of the sample exhibit distortions in their optical im- 
ages, signs of ongoing or recent interaction, and that 3.1% 
possibly have a companion. For z < 0.4, we have found that 
28.6% of the CJF sources have eith er distortions or a com- 
panion , in excellent agre ement with iHutchings & Campbelll 
(fT98l and lSuracd mm . 

- the above listed conclusion, combined with the unbiased se- 
lection criteria of the CJF sample towards transitory systems, 
imply that merger events are indeed important to the evolu- 
tion of active galaxies 

- multiwavelength variability appears to be linked to sys- 
tems undergoing a transition from a thermal to non-thermal 
emission-dominated spectrum. This may be connected to the 
presence of a bound binary black hole in the nucleus of these 
systems. 

- following the categorization of ' Lotz et al.l (l2008l) . we have 
selected the best candidate sources from our sample for dif- 
ferent evolutionary phases. We have also selected sources as 
the best candidates in our sample for BBH systems. 

More in-depth optical and infrared observations are needed, 
to unambiguously distinguish between different phenomena. 
High quality infrared data (mid and far) of our objects would 
provide us with the means to pinpoint the true position of each 



object on the assumed evolution track of galaxies. Alternatively, 
high resolution observations of those sources with multiwave- 
length variability would help us to model the possible existence 
of a binary black hole in these systems. 

Acknowledgements. M. Karouzos was supported for this research through 
a stipend from the International Max Planck Research School (IMPRS) 
for Astronomy and Astrophysics. M.K. wants to thank Jens Zuther, Mar 
Mezcua, and Chin Shin Chang for insightful discussions. The authors also 
want to acknowledge the very helpful comments from an anonymous ref- 
eree that substantially improved this paper. This research has made use of 
the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet 
Propulsion Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. This research has made use 
of NASA's Astrophysics Data Syste m Bibliographic Service s. This research has 
also made use of the VizieR service JOchsenbein et alj200(lh . 



Appendix A: Notes on individual sources 

Here we present additional information about individual sources 
from our subsample in Table |9] For references, we refer to the 
respective tables. 

Pre-merger: 



We classify as pre-merger sources systems that have a compan- 
ion but do not have distorted morphologies: 
0954+658 (BL Lac, z - 0.368): Besides its companion (at a 
projected separation 474 kpc), this source exhibits periodicity 
in the optical, and variability in y-rays. It has a bent jet on both 
VLA and VLBA scales. This is a misaligned source. The nature 
of this object is ambiguous, as it may also be classified as being 
in the post-merger phase. 

2021+614 (OW 637) (quasar, z = 0.227): Has a companion at 
a projected distance of 36 kpc. Its VLBA jet exhibits peculiar 
morphology, while its VLA jet is unresolved. 
2352+495 (OZ 488) (radio galaxy, z = 0.237): Has a companion 
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at a separation of 47 kpc. Its VLBA jet is bent, and VLA jet is 
unresolved. 



First pass: 

We classify as first pass sources systems that have distorted 
morphologies but neither a visible companion nor signs of 
ongoing starburst activity: 

0018+729 (radio galaxy, z = 0.821): This object appears to be 
distorted in the optical. On VLA scales, there is no jet resolved, 
while on VLBA scales the jet appears to be bent. 
0108+388 (OC 314) (radio galaxy, z = 0.669): Source with 
disturbed morphology in the optical. Exhibits a straight jet on 
both VLBA and VLA scales. 

0804+499 (OJ 508) (quasar, z = 1.432): Distorted morphology 
in the NIR. This source is variable in the radio. Exhibits bent 
VLBA jet and diffuse emission on VLA scales. 
1146+596 (NGC 3894) (LINER, z = 0.011): Disturbed mor- 
phology in the optical. This source is positioned in the thermal 
regime on the NIR color-color diagram (see Fig. [3]), confirming 
its transitionary phase. The above indicates that the system has 
possibly undergone, or is undergoing, a merger. The radio jet 
appears straight. A VLA jet is unresolved. 

Max separation: 



In this phase, we classify sources that show distorted morpholo- 
gies along with the existence of a companion, but no ongoing 
starburst activity. Distinguishing between pre-merger, first-pass, 
and max separation phases can be ambiguous: 
0402+379 (4C 37.11) (radio galaxy, z = 0.055): This object is 
both disturbed (in the optical) and has a close companion galaxy 
(projected distance 29 kpc). Its VLBI jet is bent. Diffuse radio 
emission is seen on VLA scales. Variable in the radio. It is posi- 
tioned in the thermal regime of the color-color plot. This object 
can also be classified as a transitory object between the max sep- 
aration and merger phase. 

0710+439 (OI 417) (Seyfert, z = 0.518): Disturbed optical mor- 
phology source with companion (projected distance 53 kpc). 
Bent VLBA jet. 

0831+557 (4C 55.16) (LINER, z = 0.24): Possibly optically dis- 
turbed source with a companion (at a separation of 81 kpc). On 
VLBA scales, the jet appears bent, while on VLA scales the jet 
and counter-jet appear misaligned. 

1031+567 (Seyfert 2, z - 0.46): Disturbed morphology in the 
NIR. A companion was found in the optical (projected separa- 
tion of 65 kpc). Shows straight VLBA jets. Unresolved VLA jet. 
1823+568 (4C 56.27) (BL Lac, z = 0.664): Source exhibiting 
non-relaxed isophotes in the NIR, indirect sign of past merger 
activity. A possible close companion is also observed (projected 
distance of 46 kpc, no redshift information). The source is vari- 
able in the radio and was claimed to show almost periodicity. 
On VLBA scales, the jet appears fairly straight, while on VLA 
scales the jet appears to be misaligned with the inner jet. 
1843+356 (radio galaxy, z - 0.764): Disturbed morphology 
source (optical) with companion (projected distance 129 kpc). 
Straight VLBA jet, unresolved on VLA scales. 
1946+708 (radio galaxy, z = 0.101): This object appears to be 
disturbed in the optical. A possible companion galaxy (no red- 
shift information available) is also observed at a projected dis- 
tance of ~ 67 kpc, thus we classify it as a max separation ob- 
ject. Its VLBA jet is highly bent. In the color-color diagram, it 
is positioned in the non-thermal regime. This may be caused by 
a misclassification. That the source is a merger remnant is the 
most plausible alternative. 



Merger: 

In the merger phase, we classify sources that exhibit distorted 
morphologies and also ongoing starburst activity: 
0248+430 (quasar, z = 1.31): This so urce exhibits signs of star- 
burst activity (e.g.. lDudik et al.|[2005h . while at the same time it 
is very luminous in both the optical and the X-rays. At a redshift 
z = 1.31, the source is not optically resolved. 
0316+413 (3C 84) (radio galaxy, z = 0.018): Disturbed mor- 
phology source (optical) with companion (projected distance 
600 kpc). This source is positioned in the thermal region of the 
color-color plot and shows ongoing starburst activity. Both the 
VLBA and VLA jet of this source appear bent. This is a primary 
candidate for a merging system. It is periodically variable in the 
radio, optical, and X-rays. 

1254+571 (Mrk 231) (quasar, z = 0.042): This object is an 
ULIRG and is positioned in the transitory region of the color- 
color diagram. It exhibits starburst activity (see Sect.|8]for more 
details). There is no jet resolved on either VLBA or VLA scales. 
Mrk 231 fits well in to the evolutionary scenario proposed by 
ISanders et al.l (1198 8') occupying a transitory phase between a 
merging system and an active galaxy. Combined with its vari- 
ability (radio, optical, and X-rays), we argue that this is a late 
merger object, moving into the post-merger phase. 
1652+398 (Mrk 501) (BL Lac, z = 0.034): This source is al- 
most periodically variable at all wavelengths checked. It exhibits 
a disturbed morphology (optical) and is positioned in the thermal 
region of the color-color diagram. This system appears to be un- 
dergoing a series of mergers, supported by it being a member 
of a cluster. Its VLBI jet is bent, with diffuse emission on VLA 
scales. 

1807+698 (3C 371) (BL Lac, z = 0.051): Optically disturbed 
source with a companion (projected separation 81 kpc). It is 
almost periodical in the optical. Variable in the infrared. The 
VLBA jet of this object appears straight. The VLA jet however 
appears bent and difiiise emission is present. 

BBH / Post-merger candidates: 



Here we select sources that exhibit variability in at least three 
wavelength regimes (almost periodic in at least one wavelength). 
Correlated variabilities are one of the strongest pieces of evi- 
dence for a BBH. These systems may also be identified as post- 
merger systems, as we expect these systems to contain a hard 
BBH system in their nuclei. 

0219+428 (3C 66A) (BL Lac, z = 0.444): This object is posi- 
tioned in the transitory region of the color-color diagram, thus 
is classified as an early post-merger system. Both the VLBI and 
VLA jet appear straight (with the VLBA jet exhibiting one sharp 
bend). It manifests almost periodic variability in the optical, and 
variability at radio and X-ray wavelengths. 
0716+714 (BL Lac, unknown z): This source is variable in 
all wavelength regimes observed, with correlated variability on 
longer timescales in both the radio and the optical. It is a mis- 
aligned source with a bent jet on both VLBA and VLA scales. 
0814+425 (OJ 425) (BL Lac, z = 0.245): Almost periodical 
variability in the radio and optical. Variable in the infrared. This 
system is positioned in the transitory region of the color-color 
diagram. Its VLBA jet appears bent. 

0923+392 (4C 39.25) (quasar, z = 0.699): This source exhibits 
variability at different wavelengths (almost periodicity in the ra- 
dio). Its VLBI jet is bent. No jet is resolved on VLA scales. It is 
positioned in the non-thermal region of the color-color diagram. 
1101+384 (Mrk 421) (BL Lac, z = 0.031): VLBA jet appears to 
be straight. The VLA jet appears bent. It is almost periodically 
variable in the radio, optical, and X-rays. It is variable in the in- 
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frared. 

1418+546 (OQ 530) (BL Lac, z = 0.151): This object has a 
disturbed morphology. That it is positioned in the transitionary 
region of the color-color diagram and it exhibits periodicities 
in different wavelengths are both indicative of a post-merger 
galaxy. Its VLBI jet does not exhibit bends. It shows almost pe- 
riodical variability in the radio. It is also variable in the infrared 
and optical. 

1633+382 (4C 38.41) (quasar, z = 1.807): Misaligned source 
with bent jet on both VLBA and VLA scales. Almost periodic in 
radio. Variable in optical and y-rays. 

1641+399 (3C 345) (quasar, z = 0.595): This source appears 
to have a close companion (at a projected separation of 47 kpc) 
and is also member of a cluster, implying the possibility of mul- 
tiple mergers. This source has already been modeled as BBH 
(iLobanov & Roland l2005h . The bent jet agrees well with the 
presence of a BBH system. Almost periodical in radio, optical, 
and X-rays. Variable also in the infrared. 

2200+420 (BL Lac) (BL Lac, z = 0.069): This system is posi- 
tioned in the non-thermal regime of the color-color diagram and 
appears almost periodical in radio and optical. It is also variable 
in infrared and y-rays. 
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